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We report a detailed and systematic study of wave propagation through a stochastic absorbing 
random medium. Stochastic absorption is modeled by introducing an attenuation constant per unit 
length a in the free propagation region of the one-dimensional disordered chain of delta function 
scatterers. The average value of the logarithm of transmission coefficient decreases linearly with 
the length of the sample. The localization length is given by £ = (, m £.a/((iw + £ a ), where £ w and 
£ a are the localization lengths in the presence of only disorder and of only absorption respectively. 
Absorption does not introduce any additional reflection in the limit of large a, i.e., reflection shows 
a monotonic decrease with a and tends to zero in the limit of a — > oo, in contrast to the behavior 
observed in case of coherent absorption. The stationary distribution of reflection coefficient agrees 
well with the analytical results obtained within random phase approximation (RPA) in a larger 
parameter space. We also emphasize the major differences between the results of stochastic and 
coherent absorption. 

PACS Numbers: 42.25. Bs, 71.55.Jv, 72.15.Rn, 05.40. +j 



I. INTRODUCTION 

Wave propagation in an active random medium has at- 
tracted much attention during the past decade. Recently, 
many experiments have reported lasing action of light in 
optically active strongly scattering media Q . These sys- 
tems exhibit interesting physical properties due to the 
combined effects of static disorder-induced multiple scat- 
tering and of coherent amplification/absorption [p|-^5[. 
In the extensively studied case of an electron motion 
in a random medium it is well established that quan- 
tum interference effects arising from a serial disorder in 
one-dimensional systems lead to Anderson localization 
|^6|,0. Studies on different types of wave propagation 
such as quantum electron transport in disordered conduc- 
tors and light propagation in random dielectric media or 
sound propagation in inhomogenous elastic media, etc. 
complement each inspite of the fact that treatment is 
quantum and classical respectively. These qualitatively 
different types of waves in an appropriate limit follow 
the same mathematical equation, namely the Helmholtz 
equation. It is basically the wave character leading to in- 
terference and diffraction which is the common operative 
feature. Light can be absorbed or amplified retaining the 
phase coherence. In most of the theoretical studies am- 
plification or absorption is modeled phenomenologically 
by introducing an imaginary potential (optical potential) 



in the Hamiltonian. In the case of light (electro-magnetic 
waves) this corresponds to a medium with a complex di- 
electric constant. Several interesting effects have been 
predicted which include statistics of super-reflection and 
transmission [^Hl3f and the dual symmetry between ab- 
sorption and amplification Media thus modeled 
are referred to as coherently absorbing or amplifying. 

In the case of electron transport, inelastic scattering 
(due to phonons) leads to loss of phase memory of the 
wave function. Thus the motion of electrons becomes 
phase incoherent and sample to sample fluctuations be- 
come self-averaging in the high temperature limit leading 
to a classical behavior. There has been much interest in 
the effect of inelastic scattering on the coherent tunnel- 
ing through potential barriers. To allow for the possi- 
bility of inelastic decay on the otherwise coherent tun- 
neling through potential barriers, several studies invoke 
absorption 



To study the above phenomenon, one 
resorts to the optical potential models (coherent absorp- 
tion models). 

In the optical potential model the potential is made 
complex V(x) — V r (x) — iVi. The Hamiltonian becomes 
non-Hermitian resulting in absorption or amplification of 
probability current depending on the sign of Vi . The pres- 
ence of imaginary potential (absorption/amplification) 
leads to many counter- intuitive features. In the scat- 
tering case, in the vicinity of the absorber, the particle 
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experiences a mismatch in the potential (being complex) 
and therefore it tries to avoid this region by enhanced 
back reflection. Imaginary potential plays a dual role 
of an absorber and a reflector |j|,|2(]]. In other words, in 
such models absorption without reflection is not possible. 
Naively one expects the absorption to increase monotoni- 
cally as a function of Vi. However, the observed behavior 
is non-monotonic |3|,po| . At first absorption increases and 
after exhibiting a maximum decreases to zero as Vi — ► oo. 
The absorber, in this limit acts as a perfect reflector. 
During each scattering event an electron picks up an ad- 
ditional scattering phase shift due to Vi which along with 
multiple interference leads to additional coherence or res- 
onances in the system Q. Thus, due to the presence of 
imaginary potentials, we have additional reflection and 
resonances in the system. In the presence of coherent 
absorption and quenched disorder, the stationary distri- 
bution for reflection coefficient has been calculated [Q. 
This has been done within random phase approximation 
(RPA) using the invariant imbedding method p2fl . The 
stationary distribution is given by: 



Ps(r) = 
= 



D\exp(\D\)exp\ 
(1-r) 2 



/ \D\ \ 

1-r for r < 1 (1) 
for r > 1. 



Here D is proportional to Vi/W, W being the strength of 
disorder. Notice that the distribution has a single peak 
which shifts towards r — with increasing absorption 
strength Vi- However, the exact distribution obtained 
numerically for strong disorder and strong absorption 
shows significant qualitative departure from this analyti- 
cal distribution j|[l3) . For sufficiently strong absorption, 
the numerically obtained stationary distribution shows a 
double peak structure. In the limit Vi —> oo the distribu- 
tion becomes a delta function at r = 1. This corresponds 
to the limit where the absorber acts as a perfect reflector. 

To this end we would like to develop a model where 
absorption does not lead to concomitant reflection and 
additional resonances as discussed above. Recently, such 
a stochastic absorption model was develo ped by Prad- 
han based on the work of Biittiker p5| , f26| . In his 

treatment several absorptive side-channels are added to 
the purely elastic channels of interest. A particle that 
once enters the absorbing or the side-channel never re- 
turns back and is physically lost. He has obtained the 
Langevin equation for the reflection amplitude R(L) for 
a random medium of length L by enlarging the S-matrix 
to include side-channels. In continuum limit the equation 
for R{L) is ||||: 

dR 

— = —aR(L) + 2ikR(L) + ikV(L)[l + R(L)} 2 , (2) 

Oiij 

where a is the absorption parameter and V(L) is the 
random potential representing the static disorder. Inter- 
estingly, within the random phase approximation (RPA) , 
the stationary probability distribution for the reflection 



coefficient P s (r) ( for L — > oo ) is again given J23 24 
by Eq.[l| In our present work we develop another simple 
model for absorption which can be readily used to study 
the case of amplifying medium as well. The medium 
comprises of random strength delta function scatterers 
at regular spatial intervals a. To model absorption (leak- 
ing out) of electrons, an attenuation constant per unit 
length a is introduced. Every time the electron traverses 
the free region between the delta scatterers, we insert 
a factor exp(-aa) in the free propagator following Ref. 
p7[ . We find that this method of modeling absorption 
does not lead to additional reflection and resonances as 
in the case of optical potential models. We obtain the 
localization length and study the statistics of reflection 
and transmission coefficients. The stationary distribu- 
tion of reflection coefficient agrees with Eq.jl in a larger 
parameter space. Following earlier method [23], the con- 
tinuum limit of our model leads to the same Langevin 
equation (Eq. ||) for R(L) where a is replaced by 2a. 
Naturally, agreement of our result with Eqjl] follows. In 
Sec.[n]we give the details of our model and the numerical 
procedure. The section after that is devoted to results 
and discussion. 



II. THE MODEL 

We carry out calculations on the wave propagation 
in an absorbing medium characterized by an attenua- 
tion constant a and interspersed by a chain of uniformly 
spaced independent delta-function scatterers of random 
strengths. The i th delta-function scattering center is de- 
scribed by a transfer matrix |2q] 



M; 



1 — iqi/2k 
iqi/2k 



-iqi/2k \ 
1 + iqi/2k J 



where qi is the strength of the i th delta-function. The 
q,'s are uniformly distributed over the range —W/2 < 
qi < W/2, i.e., P{q t ) = 1/W. Here W is the disorder 
strength. We set units of ft and 2m to be unity. The 
energy of the incident wave is E = k 2 . For further anal- 
ysis, W and a are scaled with respect to a and are made 
dimensionless. Propagation of the wave in-between two 
consecutive delta- function scatterers separated by a unit 
spacing (a — I) can be described by the matrix 



X 



nik— a 









,-ifeH 



The total transfer matrix for the L-site system is con- 
structed by repeated application of Mi and X [ pSfl : 

M = M L X....XM 2 XMi. 

From M the reflection and transmission amplitudes are 
calculated using 

M(2, 1) 



R = — 



M{2,2) 



2 



and 

detM 

T M(2, 2) ■ 

The reflection and transmission coefficients are r = |i?| 2 
and t = \T\ 2 respectively and the absorption is given by 
er = 1 — r — t. Thus, due to absorption the total flux is 
not conserved and we have r + t ^ 1 . 

III. RESULTS AND DISCUSSION 

In our studies we consider at least 10,000 realizations 
for calculating various distributions and averages. In the 
case of stationary distributions, the length of the sam- 
ples considered were about 5 to 10 times the localization 
length. We also verified that the corresponding distribu- 
tions or averages do not evolve any further with increas- 
ing sample length L. All results are shown for incident 
energy E = k 2 = 1.0 unless specified otherwise. 

We first consider the behavior of (Int). The angu- 
lar bracket denotes the ensemble average. In Fig|l] we 
plot {Int) as a function of length L for ordered absorp- 
tive medium ( W = 0.0, a — 0.05 ), ensemble av- 
eraged disordered non-absorptive medium ( W = 1.0, 
a = 0.0 ) and disordered absorptive medium ( W = 1.0, 
a = 0.05,0.1,0.15 ). In all the cases transmission decays 
exponentially with the length. The absorption-induced 
length scale £ in random medium associated with the de- 
cay of transmission coefficient is always less than both 
£ a and £ w . The localization length for disordered nonab- 
sorptive medium scales as j29| £ w — 96k 2 /W 2 and for or- 
dered absorptive medium, as £ a = 1/a. In Fig.||we show 
the plot of l/£ versus l/d w + l/£ a obtained by changing 
a for various values of disorder strength W. We have nu- 
merically calculated l/£ for the different cases. All the 
points fall on a straight line with unit slope indicating the 
relation l/£ = l/£ w + l/£ a - Such a relation exists for 
the case of coherently absorbing and amplifying media 
H[l3|. The decay of (Int) with sample length L follows 
from the general theory of random matrices also [j30| -j32| . 
In our approach scattering properties are described in the 
framework of 2x2 transfer matrices. Total tranfer matrix 
of the medium is the product of individual transfer ma- 
trices Mj of the individual scatterers. The limit L — > oo 
corresponds to multiplication of infinite number of such 
random matrices drawn independently from the same en- 
semble. In this limit the two random eigenvalues exp(±x) 
of MM* tend to the nonrandom values exp(±L / '£) with 
£ independent of L. This follows from both the Fursten- 
berg's theorem |3(J as well as the multiplicative ergodic 
theorem [|32"f . The inverse localization length 1/xi is re- 
ferred to as the Lyapunov exponent of the random ma- 
trix product in the literature. It should be noted that for 
larg e but finite L, the x has a small gaussian fluctuation 
[ p2[ around the asymptotic value £/£. This fact has an 
important bearing on the nature of fluctuations in the 
finite size sample. 



To study the nature of fluctuations in the transmission 
coefficient, in Fig.|| we plot, on log-scale, average t ((f)), 

root-mean-squared variance t v = \J (t 2 ) — (t) 2 and root- 
mean-squared relative variance t rv — t v / (t) as a function 
of length L/£ for W — 1.0 and a = 0.01. We see from 
the figure that t v is less than (t) and t rv is less than unity 
upto L/£ « 3. But, beyond that t v becomes greater than 
(t) and t rv crosses unity. We find that in the asymp- 
totic limit ln(t rv ) is positive given by the numerical value 

0. 18 L/£. The value 0.18/£ is the generalized Lyapunov 
exponent characterizing the relative variance of transmis- 
sion coefficient. Whenever the root-mean-squared vari- 
ance of a physical quantity exceeds the average value, 

1. e. when relative fluctuations become larger than unity, 
the physical quantity is said to be non-self-averaging Q . 
For example, it has been recognized long since that the 
resistance ( which is related to the transmission coeffi- 
cient of a one-dimensional random sample ex- 
hibits large statistical fluctuations when the sample size 
exceeds the localization length in absence of absorption 
or inelastic scattering (3^j35|] . The resistance fluctuations 
over the ensemble of macroscopically identical samples 
dominate the ensemble average. As a consequence the 
relative variance of resistance grows exponentially with 
the sample length L. In our study we find that in spite of 
the presence of absorption the relative variance of trans- 
mission coefficient grows exponentially with the sample 
length as mentioned earlier. Thus the transmission co- 
efficient is a non-self-averaging quantity for samples of 
length L>^. The transmission coefficient becomes very 
sensitive to spatial realizations of the impurity config- 
urations for finite size samples. This follows from the 
non-commutative nature of the Mi transfer matrices. 

In Fig.[|we show the distribution P(t) at different sam- 
ple lengths for a — 0.01 and W = 1.0. For small lengths 
L, resonant transmission dominates and P(t) peaks at 
a large value of t. In fact for L — ► 0, P(t) — > 5(t — 1). 
As the length becomes comparable to the localization 
length L ~ £, multiple reflections start dominating. Con- 
sequently, the time spent inside the medium increases 
leading to more absorption. Thus, the peak of the distri- 
bution shifts to smaller values of t and the distribution 
broadens due to randomization by disorder. In the long 
length limit L >> £, the distribution develops a long tail 
and its peak shifts towards t = 0. The transmittancc 
shows large sample-to-sample fluctuations and becomes 
a non-self-averaging quantity. Finally, as expected, for 
L -» oo , P(t) -» 5(t). 

From the previous discussion it is clear that the trans- 
mission becomes non-self-averaging. The large fluctua- 
tions in the transmission coefficient of a non-absorbing 
random medium owe their existence to the presence of 
resonant realizations (Azbel resonances) |3(],|37]]. For a 
strongly localized one-dimensional system in the absence 
of absorption at particular energies the transmission coef- 
ficient decays exponentially as a function of length L with 
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a well-defined localization length. However, for some rare 
realizations there exists a localized state close to the cen- 
ter of the sample for which incident electron can reso- 
nantly tunnel through the sample via this localized state 
with probablility approaching unity |3^,|3(|[37j . Such rare 
realizations play an important role in determining the 
fluctuations. Therefore, it is worthwhile to investigate 
the nature of resonances and the effect of absorption on 
them. Specifically, we would like to understand if the 
presence of absorption would give rise to any new res- 
onances. It is well known from the studies in passive 
disordered media that the ensemble fluctuation and the 
fluctuations for a given sample as a function of chemical 
potential or energ y a re expected to be related by some 
sort of ergodicity ||3^j58|], i.e., the measured fluctuations 
as a function of the control parameter are identical to the 
fluctuations observable by changing the impurity config- 
urations. In Fig.|^(a) we show the plot of t versus k for 
W = 1.0 and a — at L = 100 for a given realization of 
the random potential. Figure ||(b) shows a plot of t ver- 
sus k for the same realization but with a = 0.01. By mere 
visual inspection one can see that the only effect of ab- 
sorption, apart from reducing the value of transmission, 
is to increase the width of resonance peaks for the passive 
case. Thus the presence of absorption does not introduce 
any new resonances. This can be seen from Fig.^](c) and 
(d) which emphasize Fig.^|(a) and (b) respectively by en- 
larging a narrow region between k = 1.5 to A; = 2.0. We 
do not see any new peaks in the transmission spectrum 
for absorptive case. Similar effect is observed in case of 
reflection also. 

We now turn our attention to the statistics of reflec- 
tion coefficient. In Fig.|6| we plot (Inr) as a function of 
length L for a fixed value of disorder strength W = 1.0 
and different values of absorption strength a as indicated 
in the figure. It increases with L initially and for L ^> £, 
it saturates. At any L, (Inr) \w,a < (Inr) \w,o- This is 
in contrast to the behavior observed for the case of co- 
herent absorption. As we know, in the case of coherent 
absorption, the reflection coefficient tends to unity for 
absorption strength becoming very large. In this regime, 
the predominantly reflecting nature of optical potential 
makes reflection larger than that in the corresponding 
passive case. 

In Fig.]?] we have shown P s (r) for various values of a. 
In the small a range, i.e., for a = 0.001, the distribu- 
tion has a peak at large r. As we increase a the peak 
shifts to smaller values of r. The thick line shows the fit 
obtained using the analytical expression given in Eq^. 
In the limit of large a, the distribution tends to become 
a delta function at r = 0. This is in sharp contrast to 
the behavior observed for coherent absorption ||. The 
distribution is always single peaked. For all non-zero val- 
ues of a the medium acts as an absorber only and there 
is no additional reflection due to absorption. Figure ^ 
shows a monotonic decrease of saturated value of aver- 
age reflection coefficient (r) as a function of a. The 
average absorption, defined as (a) = 1 — (r) — (t), in- 



creases monotonically with increasing a and in the limit 
of a — > oo saturates to unity in contrast to the optical 
model wherein it tends to zero. In the case of the opti- 
cal model, the absorption coefficient is a non-monotonic 
function of absorption strength Vi and for values of V* 
near the peak the stationary distribution of reflection co- 
efficient displays a double peak p||n| . In fact our model 
exhibits the properties in agreement with physical expec- 
tations of an absorbing medium, i.e., stronger the absorp- 
tion lesser are the reflection and transmission across the 
medium. 

Finally, we discuss the phase distribution. Figure ^ 
shows the stationary distribution of phase of the reflected 
wave for a fixed disorder strength W = 1.0 and various 
for values of a. For small values of disorder one generally 
expects the phase distribution to be uniform if the sys- 
tem size is around the localization length. This is seen in 
Fig|](a) for the case of weak absorption. As we increase 
a the phase distribution develops two distinct peaks - 
a feature observed for coherent absorption also. This is 
related to the fact that the localization length decreases 
with a. We would like to point out that the stationary 
distribution P s (r) is same within the RPA for the case 
of coherently as well as stochastically absorbing media. 
Unlike in the case of coherently absorbing medium, Eq.|l| 
seems to be valid in a larger parameter space for the 
stochastic absorbing medium where the RPA may not 
be valid. The parameters for validity of the RPA are 
determined by the observation of uniform phase distri- 
bution. However, beyond the RPA, P s (r) for the case of 
stochastic and coherent absorbing media are qualitatively 
distinct from each other. 

In conclusion, we have studied a new phcnomcnolog- 
ical model of stochastic absorption to understand the 
statistics of quantum transport in random systems. The 
behavior observed for transmission and reflection coeffi- 
cients is in accordance with physical expectations of an 
absorbing medium. This model can be extended to the 
case of stochastically amplifying medium. It exhibits du- 
ality between absorption and amplification which has re- 
ceived much attention recently. Results for this will be re- 
ported elsewhere . It is to be noted that the treatment 
is a phenomenological one. A better treatment based on 
first principles like density matrix involving system and 
its coupling with environment is called for. 
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FIG. 1. Average of logarithm of transmission coefficient t 
versus length L for different values of absorption strength a 
and disorder strength W as indicated in the figure. 




FIG. 2. l/£ versus l/£ w + l/£a, where £ TO is the localiza- 
tion length for non-absorptive disordered system, £ a is the 
localization length for ordered absorptive system and £ is the 
localization length for absorptive disordered system. 
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FIG. 3. Average transmission coefficient (t), variance of t 
and relative variance of t as a function of L/£ for fixed disor- 
der W = 1.0 and absorption a = 0.01. 
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FIG. 5. Transmission coefficient as a function of incident 
wavenumber k for (a),(c) disordered non- absorptive sample 
W — 1.0 of length L = 100 and (b),(d) disordered absorptive 
sample W = 1.0, a = 0.01 of length L = 100. 
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FIG. 4. Distribution of transmission coefficient t from a 
disordered absorptive system with W = 1.0 and a = 0.01 at 
different lengths L. 




FIG. 6. Average of logarithm of reflection coefficient ver- 
sus sample length for a fixed value of disorder W = 1.0 and 
different values of absorption a indicated in the figure. 
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FIG. 7. Stationary distribution of reflection coefficient for 
a fixed value of disorder strength W = 1.0 and different val- 
ues of absorption a. The thick line shows the single parame- 
ter fit of analytical expression Eqn.|l] with (a) D = 0.197 for 
a = 0.001, (b) D = 1.92 for a = 0.01 and (c) D = 20.53 for 
a = 0.1. 
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FIG. 8. Average value of reflection coefficient ((r)), trans- 
mission coefficient ((i)) and absorption ({a}) versus absorp- 
tion strength a for W = 1.0 and L/£ = 10. 
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FIG. 9. Stationary distribution of phase of reflected wave 
for fixed disorder strength W — 1.0 and different values of 
absorption a. 
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